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Abstract 



In spite of their evident logical character, particle statistics symmetries are 
not among the inherently quantum features exploited in quantum compu- 
tation. A difficulty may be that, being a constant of motion of a unitary 
evolution, a particle statistics symmetry cannot affect the course of such an 
evolution. We try to avoid this possible deadlock by introducing a generalized 
(counterfactual, blunt) formulation where this type of symmetry becomes a 
watchdog effects shaping the evolution of a unitary computation process. The 
work is an exploration. 

I. INTRODUCTION 

This study deals with an alternative form of quantum computation. In this context, it 
is useful to identify the "boundary" of the usual approach. 

Until now, quantum computation speed up^^"^ among others] ^^j^]-^ respect to the known 
classical algorithms) has taken advantage of a wide but incomplete set of inherently quantum 
features: effective reversibility, superposition, entanglement, interference and measurement. 
However, particle (fermionic, bosonic) statistics has never been exploited. 

We should further note that all quantum algorithms found so far are based on the same 
form of computation, which is time- sequential, in the sense that a reversible Boolean network 
appears in the time-diagram of the logical transformations undergone by the register. 

They are also based on a common strategy, as recently recognizedl^"^' ; for example, they 
can be seen as sophisticated applications of the Mach-Zender interferometer. 

Furtherly, this form of computation is suspected'^' to be inherently unable of solving in 
polynomial time an NP-complete problem.^ 



Solving in polynomial time (P) a problem formerly classified in NP (in particular, whose solution 
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The alternative approach explored in this study is tailored to solve a particular NP- 
complete problem'^^l and is essentially based on the exploitation of particle statistics - 
under a generalized, counterfactual and blunt interpretation thereof. 

We give up time-sequential computation, by considering, in place of the register, a re- 
versible Boolean network fully deployed in space (rather than appearing in the register 
time-diagram). Time is now orthogonal to the network lay-out (fig. 1). In this way, part of 
the network input and part of the output can both be constrained (fig. 1). Checking whether 
this network is satisfiable is a version of the NP-complete satisfiabihty (SAT) problem. 

We shall describe this in some detail; in the meantime, given the novelty of the approach, 
we give an outline of the entire work. 



Fig. 1 
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• a,b, ...,h (fig.l) labeling the network nodes, are Boolean variables. The overall network, 
or each gate, has its own input (left) and output (right) nodes- variables: a, b, f (c, g, h) 



required exponential time with all formerly known algorithms), means declassifying to P that 
problem alone. Solving in polynomial time an NP-complete problem would declassify to P all 
NP-complete and all NP problems. 
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make the input (output) of the overall network; a, b (c, d) make the input (output) of 
gate #1, etc. With time orthogonal to the network lay-out, the terms "input" and 
"output" loose part of their meaning but are kept out of habit. 

• Each gate introduces a logical constraint between its input and output variables, by 
establishing a function between them. This function is given in table form: the gate 
"truth table". For example table 1, associated with gate 1, gives the invertible XOR 
(exclusive or) function. Table 2, associated with gate 2, gives the NOT function (to 
spare notation, the two-node NOT gate will be used as a network wire - an invert- 
ing wire, without loss of generality). Tables 4, 5, 6 represent the (partial) network 
input/output constraints. 

• Table N is introduced for convenience; it does not estabhsh any constraint, its 2^ rows 
are all the possible assignment of the network Boolean variables. 

• The SAT problem is whether there is at least one Boolean assignment - one row of 
N - satisfying all the network constraints (all gates and wires and the input /output 
constraints) . 

• The above is readily translated in Hilbert space mathematics^^^l . A specific assign- 
ment of all network Boolean variables (one row of N) becomes a tensor product of 
corresponding qubit eigenstates, and a basis vector of the network Hilbert space Hn- 
The logical constraint established by each gate or wire (network element or N.E.) i can 
be seen as a projector Ai from Hn on Ti^f^. l-6f' is a "constrained" Hilbert subspace, 
spanned by those basis vectors of Tiiv which satisfy the constraint established by N.E. 
i. 

• Qubits are implemented as a distributed system of two-state identical particles. The 
particle distribution corresponding to the input/output nodes of N.E. i is subject to 
a Hamiltonian H^. Given a suitable form of ifj, the projector is physically "turned 
on" in the degenerate ground state of Hi, more generally in any network state satisfying 
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-ffj = 0. Such projectors are shown to be epiphenomena of identical particle 
statistics, under the generalized interpretation'^^'^^^. All Hi commute pairwise and the 
network Hamiltonian is if at = Yli^i- 

• The network is prepared in a specific ground state I^E' (0)) of if at 0. {'^ (0)| if at I^P (0)) = 
implies (0)| ifj |^ (0)) = for all i, consequently all Ai are "turned on" in |^ (0)). 
Thus |\E'(0)) satisfies IljAj (0)) = |\I^(0)), namely all N.E. constraints (all Ai are 
pairwise commuting). This preparation is also a linear combination of all basis 
vectors which satisfy the network input constraints, not necessarily the output con- 
straints. Since the output constraints are removed in this phase, until now computation 
has solved a simple P (polynomial) problem. 

• By suitably operating on the output qubits in order to bring them in match with 
their constraints, the network ground state, continuously submitted to the projection 
IliAi \ (t)) = |\E' (t)), is unitarily brought from |^E' (0)) to a state |\t'' (r)) satisfying 
also the output constraints (if the network is satisfiable). This requires interpreting 
a particle statistics symmetry as a continuous watchdog effect shaping the network 
evolution It should be noted that (r)) can be in superposition with an orthog- 
onal excited state |\t'" (t)), which is an error state from a computational standpoint. 
However, the amplitude of this latter state should be small. 

^The generic ground state of if at is 2" degenerate, where n is the number of input (or output) 
qubits. 

^An example of a unitary evolution shaped by a watchdog effect, is the evolution of the polarization 
of a photon going through an infinite series of polarizing filters, each rotated by an infinitesimal 
constant angle with respect to the former one. In a way, we go back to the root of quantum 
computation (computation reversibilityt^'^'^^]) and take an alternative branch, by exploring a strictly 
quantum form of reversible computation. 
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• Measuring the network after the above operations gives a Boolean assignment. If 
the network is satisfiable, this is hkely to be a network solution (whether it is a so- 
lution is checkable off line in polynomial time). If the check gives a non-solution, 
the network is likely to be not satisfiable. By repeating the overall process for a 
sufficient number of times, it can be ascertained with any confidence level whether 
the network is satisfiable. Computation has now solved the hard, NP-complete, 
SAT problem. By the way, the constrained network of fig. 1 is satisfied by (only) 

l^) = |l)a|l)Jl)c|0)Jl)e|0)/|l),|l).. 

• The possibihty that particle statistics yields a computation speed-up appears to be 
promising in principle. 

II. A COUNTERFACTUAL INTERPRETATION OF PARTICLE STATISTICS 

This work hinges on a counterfactual interpretation of particle statistics^^^'^^'^^'. A parti- 
cle statistics symmetry is usually considered to be an initial condition which is conserved as 
a constant of motion. This notion can be made more general by applying a counterfactual 
logic: the state of the system might violate the symmetry, but it cannot be so because of 
a special watchdog effect which would immediately project such a state on the symmetric 
subspace. This way of reasoning yields to effective consequences^, diverging from the con- 
ventional application of quantum mechanics. Apropos of this, it is worth citing the following 
statement by Roger Penrose: "What is particularly curious about quantum theory is that 
there can be actual physical effects arising from what philosophers refer to as counter f actuals 
— that is things that might have happened, although did not in fact happen" t^^l .|^ 

^besides the fact that particle statistics symmetries should no more be put among the initial 
conditions. 

^The relevance of counterfactual reasoning in sequential computation has been highlighted by 
Richard Jozsa^^'^^. 
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By way of exemplification, let us consider a state which we assume to be symmetrical 
(for particle permutation) because of particle statistics: 



*(0))=cos^^|0), 10)2 + sin^cos^(|0)Jl)2 + |l)i|0y+sin2^|l), 11)2; (1) 



1 and 2 label two identical two-state particles (which makes this an idealized situation) . The 

above amplitudes arc used for convenience with no significant restriction. Of course |^ (0)) 
is symmetrical under the particle permutation P12 iff it satisfies the equation 5*12 |^ (0)) = 
1^(0)), where S'12 = | (1 + P12) is the symmetrization projector. Let Hu be the Hilbert 
space of the two particles and T-Ci2 be the symmetric subspace: 



In order to introduce the counterf actual picture, we assume that |^ (i)) might violate 
the symmetry. This is interpreted by taking, for all t, a |^ (t)) which is a free normalized 
vector of Ti.i2- Namely, given: 



the aij are free complex variables independent of each other up to normalization for any 
time t. In other words, if 7^ ^2, |^ (^1)) and I^P (^2)) are two independent free normalized 
vectors of Hu- 

The logic is that, at any time t — t-, |^ (t_)) might violate the symmetry: it is thus a 
free vector of lH-i2- Then |\E' (t_)) should be projected on li-12. This gives as a result the 
vector of 1-612 closest to |^ (i_)) . 

Let 1^ {t+)) be a free normahzed vector of 7Yi2 independent from |^ {tJ)) {t+ > t-). The 
above projection can be obtained by submitting |\1' (t_|_)) to the following conditions: 






ije{o,i} 



ije{o,i} 



for all t: 



i) S,2\^{t+))^\^{t+)), 



ii) |(^ (i+) I ^ (^-))| must be maximum. 
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In this picture, (t-)) is the vector "immediately before projection" while |\E' subject 
to conditions (i) and (ii), is the vector "immediately after projection". This counterfactual 
logic will become clearer when applied to an evolution of the above triplet state. 

Continuous projection can be interpreted as a continuous state vector reduction induced 
by particle statistics. To show this, it is convenient to rewrite (2) as follows: 

1^ (t-)) = «io"^ |0)i 10)2 + I {aV + c^V) (|0)i 11)2 + 10)2) + 
I {c^V - c^V) (|0)i 11)2 - 10)2) + 11)2 • 

This vector undergoes a partial ^ state vector reduction on the subspace ■ This cancels 
the amplitude of the channel (|0)^ |1)2 — 11)^ 10)2) and renormalizes the amplitudes of the 
other channels. This can be seen as a special form of interference (from now on: "interfer- 
ence" ) , destructive in amplitude cancelation and constructive in renormalization. Of course 
this form of state vector reduction is an interpretation of particle statistics and is not due 
to an interaction with an external apparatus. 

t+ might be as close to t_. as desired, since state vector reduction does not involve the 
notion of elapsed time. However, we will take t_ = t — dt and t+ = t, in view of applying 
(i) and (ii) to an evolution described by \^ {t)). The former conditions are thus rewritten: 

for all t: 

i) 5i2|vI/(t)) = |M/(t)), 

ii) Max |(^(t) I ^{t-dt))\, 

where \^ {t)) is a free normalized vector of 'Hi2- 

Now we go back to our system in a triplet state. Duly constrained, 1^1' (t)) will describe 
an evolution of that state. Conditions (i) and (ii) will be applied together with an operation 



f\ ' is] 

It is essential that 'H12 has dimension higher than one. 
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performed on the system, say an identical rotation of the state of both particles of an angle 
ip (t) (where ip (t) is some function of time), starting from preparation (1). 

Let us anticipate the result, by following the conventional approach. The above oper- 
ation, symmetrical for particles permutation, commutes with 5'i2, therefore symmetry is a 
constant of motion of the evolution; namely, it does not apply any constraint to it. The two 
particles can thus be rotated independently of each other. This yields of course: 

I* it)) = cos' l^ + p it)] \0), \0), + sin i^ + p it)] cos + p {t)] (|0), |1)2 + |l)i 10)2) (3) 
+ sin2 [^ + p>{t)]\l),\l). 

Now we want to obtain the same result by ignoring the fact that symmetry is conserved as 
a constant of motion, and by resorting to the projection interpretation of particle statistics, 
that is to conditions (i) and (ii). We need to add a further condition to specify that the 
state of each particle has been rotated by cp (t) with respect to the initial state (1). 

Doing this requires some attentions. Now the two particles cannot be considered inde- 
pendently of each other, since their overall state is submitted to the continuous projection, 
or watchdog effect, (i) and (ii). We must utihze the particle density matrices: 

A(i)=rr3_,[l*W)(*WI], 

where Tr^^i means partial trace over 3 — i, and i = 1, 2 is the particle label (if we use the 
method of random phases'^^l, (t)) does not need to be a pure state - anyhow it will turn 
out to be that). Furthermore, the coherence elements of each density matrix - as overall 
entanglement - can be affected by the watchdog effect (in fact they will be determined by 
it). We only know that the diagonal of each density matrix must show a ip (t) rotation. In 
conclusion the further condition (iii) is: 

diagp^ (t) = diag {Tr^ [|* {t)) (t)|]} = cos^ + p {t)] (Oj^ + sin^ + p {t)] (Ij^ . 
diagp2 (t) = diag {Tn [|^ (t)) (t) |]} = cos^ + p it)] (Oj^ + sin^ + p {t)] . 
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It is readily seen that the simultaneous apphcation of conditions (i) and (iii) yields: 
\^ (t)) = cos^ + ^ (t)] |0)i 10)2 + e^'^ sin [{} + ^ (t)] cos + ^ (t)] {\0), + \0),) 

where the phases 61 and 62 are still unconstrained; condition (ii), given the initial state (1), 
sets 61 = 62 = and keeps this result frozen throughout the evolution, as readily checked. 
This yields the desired evolution (3). 

As readily seen, by removing either one of the two conditions (iii), evolution (3) is still 
obtained - the two rotations are redundant with respect to one another. 

Therefore it is perfectly legitimate to say that the rotation of the state of only one particle, 
either one in an indistinguishable way, drags an identical rotation of the state of the other. 

It should be noted that applying condition (iii) to either part of the system does not mean 
applying the identity to the other partQ. In the current formulation, the transformation of 
this latter part must be disregarded, since it is already determined by the transformation of 
the former part and the watchdog effect. In a sense, either part, being indistinguishable, 
represents the whole. This might seem a peculiar way of applying quantum mechanics but 
it is tautologically true here. 

In this idealized picture (the two-state particles), particle statistics symmetry can be 
seen as an interaction-free constraint (working by way of "interference") between the states 
of the two particles, before being a constant of motion. Of course symmetry is conserved by 
the propagation (3), but this is due to the fact that the propagation is shaped in a symmetric 
way by the watchdog effect described by (i) and (ii). 

This work explores the possibility of porting the above, generalized, interpretation of 
particle statistics to the case of identical particles hosted by distinguishable sites, provided 



^Which would mean: (a) applying an infinitesimal rotation to one particle and the identity to the 
other, (b) symmetrizing the result and (c) repeating the cycle. This would bring in a Zeno effect 
freezing the evolution in its initial state (1) - see ref. [12]. 
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that the whole system is in the ground state of a suitable Hamiltonian. Section III gives a 
purely mathematical porting. The physical model will be the subject of Section IV. 



III. THE NOT GATE AS A PROJECTOR 

We consider a simple NOT gate (fig. 1, gate and table #2). When used as a network 
(inverting) wire it will be called "link" . 

Mathematically, the previous model is readily ported to a pair of qubits r and s [rs would 
run over de, etc. in a generic network] , of Hilbert space 

Hrs = span{\0)^ |0), , |0), , |0), , |1) J. 

These two qubits will be implemented in a way that makes them distinguishable (Section 
IV). 

The constrained subspace, spanned by the basis vectors of Hrs which satisfy the NOT 
gate, is: 

= span {|0)J1),,|1) JO) J. 
The projector Ars from Hrs on H^} is defined by the following equations: 

^r.|0)Jl), =|0)J1),, A,,|1)J0),= |1)J0),, 

A.|0)J0), =0, A,|1)J1), = 0. 

Finding a physical model of Ars will be the subject of the next section. 

The two qubits should be prepared in an initial state already belonging to H^]: 

|*(0))=cos^|0)Jl), + sin^|l)jO),. (4) 

The evolution we are looking for, supposedly obtained by operating only on qubit r and 
disregarding qubit s (see Section IV) , is the solution of the following system of simultaneous 
equations: 
For all t: 
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i) A.|*W> = |*W), 



ii) Max \{^{t) \^{t-dt))\, 



iii) diagpr (t) = diag {Tr, [|* {t)) {t)\]} = cos^ + ^ (t)] |0)^ (0|^ W + if (t)] , 



where |^ (t)) is a free normalized vector of H.^ 

2 

= 1. 



= Dje{o,i} "S? N)r U). ' 



Q;. • 



Zli,ie{o,i} 

As readily seen, conditions (i) and (iii) yield: 



I* (t)) = cos [^ + ^ it)] |0), + e^' sin + ^ (t)] |0), , 

where 6 is an unconstrained phase; condition (ii), given the link initial state (4), sets 5 = 
and keeps it frozen throughout the evolution. This yields to the unitary evolution: 



I* (t)) = cos 1^ + ^ it)] |0), |1), + sin [^ + <p it)] |1), |0), 



(5) 



We must assume that condition (iii) represents some physical operation performed only 
on qubit r, under conditions (i) and (ii) - exactly like in Section II. The plausibility of this 
assumption will be discussed in Section IV. This would mean that a transformation operated 
on qubit r drags an identical transformation of qubit s. Mathematically, this is true: 



Ps (t) = Tr^ [|* (t)) m = sin^ + (t)] 



with respect to condition (iii), and 1 are interchanged because qubit s is the NOT of qubit 



r. 



Interestingly, if all Til^J basis vectors occur with amplitudes different from zero in the 
initial state (4), namely if 7^ 0, |, condition (i) is redundant with respect to condition (ii). 
In this case, condition (ii) alone implies ccqq = af^ — 0, and this already satisfies condition 

On the contrary, condition (i) is not redundant if = or |. For example, if = 0, i.e. 



1^(0)) = 10)^11)5, condition (ii) implies a^Q — and 



a 



cos^ [& + (p{t)], as needed. 
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a 



(t) 
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+ 



a 



it) 
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sm"^ [d + if (t)]. Thus, 



while afj and afi are only subject to the constraint 
disregarding condition (i) would allow for the existence of the forbidden state afl 

It is worth noting that the same evolution (5) can be obtained by applying the unitary 
operator Qrs [^p {t)] to the overall state |\E' 



Qrs b (t)] = 



^10 0^ 

COS ip (t) sin ip (t) 
— sin (t) cos ip (t) 















with |0)J1),^ 



W 



|1)J0),= 



W 




Qrs (t)] operates on the overall state in a non factorizable way, since at any time t the 
transformation undergone by qubit r (s) is conditioned by the state of qubit s (r), because 
of entanglement: 

Qrs (t)] |0), = cos^ (t) |0), + sin^ (t) |0), , 

Qrs b W] |0), = - sm^it) |0), + cosy, (t) |0), . 

There is so to speak a "hidden" conditional logic: finding this evolution — i.e. Qrs [<p {t)] 
— amounts to solving a logical problem. This becomes NP-complete when the link or the 
NOT gate belongs to a generic Boolean network. 

We have thus ascertained a peculiar fact. Our "operation on a part" [this is just the 
mathematical condition (iii) for the time being] , blind to its effect on the whole, performed 
together with continuous Ars projection, generates a unitary transformation which is, so to 
speak, wise to the whole state, to how it should be transformed without ever violating Ars 
(i.e. the link or the NOT gate). Of course Ars ends up commuting with the resulting overall 
unitary propagator, but because this is shaped by Ars- 



*the generator of Qrs [^{t)] commutes with Hrs, the interaction Hamiltonian between qubits r 
and s which will be introduced in the next Section. 
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IV. TURNING ON THE NOT GATE 



projection can be shown to be an epiphenomenon of particle (fermionic or bosonic) 
statistics "turned on" in a special physical situation. In the following, we will adopt fermionic 
statistics. 

We submit a couple of identical fermions 1 and 2 to a suitable Hamiltonian. Each 
fermion must have two compatible, binary degrees of freedom x and ^- Just for the sake 
of visualization (things should remain more abstract), we can think that each fermion is 
a spin 1/2 particle which can occupy one of either two sites of a spatial lattice, x can be 
the particle spin (say that x = 0, 1 correspond to spin =down, up) and X — r,s the label 
of the site occupied by the particle. For example, \0)^ |1)2 |r)^ |s)2 reads: particle 1 spin = 
down (0), particle 2 spin = up (1), particle 1 site = r, particle 2 site = s. There are 16 
combinations hke this, which make up the basis of the Hilbert space Tiu- However, there 
are only six antisymmetrical combinations (not violating fermion statistics) which make up 
the basis of the antisymmetrical subspace . 

These basis vectors are represented in first and second quantization and, when there is 
exactly one particle per site, in qubit notation {x and A stand respectively for the qubit 
eigenvalue and label), |0) is the vacuum vector: 



k) = 75 (I0)i 11)2 - 10)2) k)i \r)2 = «l «l |0) , 
1^) = ^ (|0)i 11)2 - 10)2) 1^)2 = al al |0) ; 
|c) = 71 |0)i 10)2 (|r)i 1^)2 - \s), \r),) = al al |0) = |0), |0), , 

\d) = 75 11)2 1=^)2 - k)2) = «l 4. |o) = , 

|e) = Hlo)i 11)2 + 10)2) k)2 - k)i 1^)2) = 

^ (al al + al al) |0) = ^ (|0), |1), + |1), |0) J . 

I/) = |(|0)i 11)2 - 10)2) k)2 + k)i |r)2) = 
^ (al al - al al) |0) = ^ (|0), |1), - |1), |0)J . 



a* . creates a particle of spin x in site A; creation/annihilation operators are subject to: 




Under the condition that there is exactly one 
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particle per site, they generate a qubit algebra. 
Now we introduce the Hamiltonian 

Hrs = Ea \a) {a\ + \h) {h\ + E^ \c) (c| + E^ \d) {d\ 

or, in second quantization, 

Hrs = -{Ea a\)r al^aorOir + -^6 oj^ aLoosOis + Ec al^ al^aoraos + Ed aLoirOis), 

with Ea, Eb, Ec, Ed > E discretely above 0. This leaves us with two degenerate ground 
eigenstates: 

|e) = ^ (|0), + |0) J and |/) = i= (|0), - |0)J . 

The generic ground state is thus: 

l*) = «|0)Jl), + /3|l)JO),,with + = (6) 

Of course |^) satisfies Ars |^) — |^), and belongs to H,^} (Section III), a subspace of Ti!f2 ■ 
Let Ai2 \^) — \ {1 — P12) be the usual antisymmetrization projector. Due to the above 
anticommutation relations: 

A2|0)J1), = |0)J1), andAi2|l),|0), = |l)jO),, moreover 

|0), |0), = |0), |0), and |1) Jl), = , 

without forgetting that |0)^ |0)g = |c) and = \d) are excited states. 

The link can be implemented by suitably operating on the ground state (6). We assume 
this to be initially given by eq. (4), repeated here: 

1^ (0)) = cos^9 |0)^ + sin^? |0)^ . 

The transformation of qubit r under (the equivalent of) a continuous Aj-s projection is 

described by: 
for all t: 
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i) Ai2|*W) = |*(t)), 

ii) Max \{^{t) \ ^{t-dt))\, 

iii) diag pr {t) = diag{Trs [(|* {t)) {t)\)]] = 
cos' [i} + cp (t)] |0), (0|, + sin^ + ^ (t)] , 

iv) {irAt)) = {'^{t)\Hrs\^{t))=0, 

where |^ (i)) is a free normalized vector of Hi2- 

The solution of the above equations is the desired link evolution (5), repeated here: 

I* it)) = cos i^ + ip it)] |0), + sin i^ + ip (t)] |0), . 

Conditions (i), (ii) and (iii) mean that the link state undergoes a transformation [driven 
by (iii)] under continuous state vector reduction on the antisymmetric subspace . If 
■& 7^ 0, |, namely if the preparation (4) comprises all the basis vectors of Hr'^s, condition 
(ii) alone keeps the link evolution inside Tii'^J (Section III). The link state remains ground 
and consequently the link expected energy {^rs (t)) is always zero. By excluding i? = 0, |, 
condition (iv) is a consequence of the former conditions. 

Mathematically, conditions (i) and (iv) give the constraint Ars |^ (^)) = |^ (^))- In 
conclusion the above conditions (i) through (iv) (which imply interpreting fermionic an- 
tisymmetry Ai2 as continuous projection on the antisymmetric subspace) are equivalent to 
condition (i) through (iii) of Section III. This gives in fact the evolution (5). 

Let us see how we could operate (in principle) on the state of qubit r, in order to im- 
plement the driving condition (iii). Since conditions (i) through (iv) generate (under the 
counterfactual interpretation of particle statistics) the projector Ars, operating on qubits r 
and s, the state of qubit s is redundant with respect to the state of qubit r. In particular, 
operating only on qubit r does not mean applying the identity to qubit s, it means disregard- 
ing it (we are now inside that "peculiar way of applying quantum mechanics" highlighted 
in Section II). 
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In order to operate on qubit r, we introduce the one-qubit Hamiltonian Hr — E'raJr^or of 
ground state and excited state lO)^.. Hr and Hrs commute and the total expected energy 
of the system is {t) \ {Hrs + Hr) \ ^ (t))- This is conveniently split into a link internal en- 
ergy {^rs (t)) = (t) l-ffrsi * (t)) and a qubit r internal energy (^^ (t)) = {"If (t) * (i)) = 
Er cos^ [■(9 + (/? (t)]. One can see that the latter is in one-to-one correspondence with diag 
Pr {t): the driving condition (iii) can thus be implemented by changing the internal energy 
of qubit r, provided that {^rs (t)) = (see further below). 

Qubit r is now put in interaction with a heat-bath. This could also excite the link, 
yielding {^rs {t)) = {t) \Hrs \ ^ {t)) > 0. In order to take this possibility into account, it is 
convenient to split l^' (i)) into two orthogonal components: 

\^(t))^a(t)\^'(t))+P(t)\^"{t)), 

with la {t)f + \(3 {t)f = 1, such that (t) \Hrs\ (t)) = 0, (t) \Hrs\ (t)) = {^rs (t)) > 
0, \'^' (t)) belongs to hI'^s and 1^^" (t)) belongs to the orthogonal subspace. 

The former component is hosted, so to speak, in a "good" Everett universe where 
{Crs {t)) — 0, Ars is turned on and the link works as required: the transformation per- 
formed on qubit r drags a corresponding transformation of qubit s. Prom the standpoint of 
the link internal energy, this occurs in a reversible way with no free energy dissipation (nor 
increase) . 

The latter component is hosted in a "bad" Everett universe where the link state is a 
linear combination of the excited states \0)^ |0)^ and and the link (or NOT gate) 

logical constraint is not satisfied. 

P (0) was zero in the preparation (4). Let T be the heat-bath temperature. We assume 
that ksT << E ^ Ea, Ef,, etc.. Thus (3 (t) will remain close to zero after the heat-bath has 
been turned on. Therefore, the probability |Q;(i)|^ that measurement finds the link in the 
"good" Everett universe should remain very high. 

We further assume that E^ << E is of the order of ksT. {C,r (t)) is driven to relax close to 
zero by reducing T. At a time r such that {^r (t)) = cos^ [t!} + (p (r)] ~ 0, and in the "good" 
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Everett universe, the link state has the form (r)) ~ \0)^. The time required to reach 
{^r (t)) — is uniquely determined by the relaxation process of qubit r - independently of 
the rest of the link. 

This way of driving the link evolution, whose plausibility is based on heuristics for the 
time being, will be used to drive the evolution of the entire network (Section VI). 

V. THE GENERIC GATE 

To reduce notation, we shall work with a gate of three (coexisting) qubits, the minimum 
required to have a conditional logic. The result will be clearly generalizable to any number 
of qubits. We will consider the logically irreversible XOR gate, of inputs t, u, output v, and 
thruth table: 



I 





t 


u 


V 














1 


1 


1 





1 


1 


1 






Since time is orthogonal to the network lay-out, the network evolution can be reversible 
even if gates are logically irreversible. 

The gate mathematical model will be introduced first. Let 

Htuv = span{\0), |0)„ |0), , |0), |0), , |1) Jl) J 

be the eight-dimensional Hilbert space of the gate three qubits. The constrained subspace 
(whose basis vectors map the rows of the gate truth table) is four-dimensional: 

^£ = span m, |0)„ |0), , |0), |1)„ |1), , |1), |0)„ |1), , |1), |1)„ |0) J ; 

The projector Atuv from Tituv on 7i[uv is defined by the (eight) equations: 
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Auv \Xa)t Mu IXt)^ = \Xa)t IXt)^ when XaXpXj is a row of the gate truth table, 
Auv \Xa)t \Xp)u Mv = when XaXpXj is not a row of the gate truth table (like e.g. 

tuv. 

010). 

Let us denote by (t)) the evolution of the gate state. (t)) is continuously projected 
on 7i[uv This is represented by: 

for all t : 

i) Atuv\'^{t)) = \^{t)), 

ii) Max \{^{t) \^{t-dt))\, 

where |^ (i)) is an otherwise free normalized vector of Htuv', in fact conditions (i) and 
(ii) do not yet determine |^ (t)): the "missing" driving condition (iii) will come out 
from the behaviour of the rest of the network (Section VI). 

The extension to a gate with any number of qubits is obvious. 

Now we should build the gate model Hamiltonian Htuv We need three lattice sites 
t,u,v and three identical spin | particles 1, 2, 3. The unwanted states should all be excited 
eigenstates of Htuv This has therefore the form: 

Htuv — ■■■ ~ Ei^aQ^aitOQyaotCiitO'Qv ■■■ — E^aQtOiuf^ov^Qt^^iu^^Qv ■■■ 

The first term corresponds to a state where each site is not occupied by exactly one particle. 
The second term corresponds to a state where each site is occupied by exactly one particle 
but is not a basis vector of H^uv^ ^tc. 

Conditions (i), (ii) and (iv) of Section IV become now: 

for all t: 

i) Ai23 1^ (t)) = 1^ (t)), where A123 = ^ (1 - P12 - P13 - P23 + P13P12 + P23A2) is the 
antisymmetrization projector, 
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ii) Max|(*(i) 

iv) (Ctuvit)) ^ {t)\Htuv\^ {t)) ^0, 

where (t)) is an otherwise free normahzed vector of 7ii23 (similar to of Section 
IV). 

The gate initial state should be a linear combination of all H^^y basis vectors. Consequently, 
because of (ii), the gate state will remain in Hfuv- condition (iv) is thus a consequence of 
(i) and (ii). (i) and (iv) are equivalent to Atuv |^ (^)) = 1^ (^))- ^tuv is thus "turned on" 
when {^tuv (^)) = 0. 

VI. THE NETWORK 

We must put together the stand-alone network elements of the previous Sections. The 
network lattice is the union of the lattices of the individual N.E. Of course not all sites 
are disjoint: for example a gate and a connected hnk share one lattice site. All N.E. 
Hamiltonians Hi (where i is a multiple index running over de, abed, etc. - fig. 1) are 
diagonal in Hn basis thus pairwise commuting: Vjj : [Hi, Hj] — 0. 

The overall network Hamiltonian is: ifjv = ^^'^ network expected energy 

(Ctv (t)) is the sum of the expected energies of all N.E.: {^n (t)) = (^)l 1^ (^)) > 

where |^ {t)) is the network state (i.e. a linear combination of the network basis vectors - 
each a tensor product of all qubit eigenstates) . 

Therefore, if |\E' (t)) is a nework ground state (this is highly degenerate, see Section I), 
then (^jv (^)) = and consequently {t) \ Hi\^ (t)) — for all i. Fermion statistics and 
these latter conditions work as conditions (i) and (iv) of Section V: all Ai projectors are 
"turned on" when (^^v (t)) = 0. 

Of course Vj^ : [74j,>lj] — 0. An — HiAi projects Hn on the constrained subspace 
spanned by the network basis vectors which satisfy all N.E. (Hn' is the intersection of all 
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To sum up, any transformation of the network ground state (t)) satisfies the following 
conditions: 
for all t : 

i) A^|vl/(t)) = |v^(t)), 

ii) Max \{^{t) \ ^{t-dt))\, 

where l"^ (t)) is an otherwise free vector of Hat- 

Of course, any tensor product of the network qubit eigenstates satisfying condition (i) sat- 
isfies all the N.E., not necessarily the input and output constraints. 

The satisfiability problem requires checking whether such constraints can also be satis- 
fied. Under the assumption that the network is satisfiable together with the input-output 
constraints, we proceed as follows. 

1) We assume, with no restriction, that only one output qubit, say z, is constrained: 
Pz = \l)z (Mz- This output constraint is temporarily removed. Whereas we introduce one- 
qubit Hamiltonians in order to detain the constrained inputs in the required values. Let h 
be an input qubit constrained by = |0)^ (0|^ {ph = (Mh)- This constraint is satisfied 
in the ground state of the one-qubit Hamiltonian Hh = Eha\^aih {Hh = Ehalf^aoh), with 
Eh discretely above zero (of the same order of E). We should note that these Hamiltonians 
commute with all the other N.E. Hamiltonians. 

2) We wish to prepare the network in a state |\E' (0)) which is a linear combination of 
all T-Cn basis vectors satisfying the input constraints, all the N.E., but not necessarily the 



^We should note that the permutation of two identical particles belonging to any two different N.E. 
cannot set any further particle statistics constraint, since the two N.E. Hamiltonians commute. 
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output constraint (not implemented until now). This means satisfying the equations: 

• V /i : {0) \ Hfi]^ (0)), where h ranges over the labels of the input qubit and N.E. 
Hamiltonians; 

• y k : pk (0) — Qk \0)i^ {0\i^ + bk , with a^, bk ^ 0, where k ranges over the labels 
of the unconstrained input qubits; we should note that the coherence elements of any 
qubit density matrix are zero if, without any restriction, any two qubits are connected 
by a link (although this is not the case of fig. 1). 

The above conditions yield: 

I* (0)) = cos Ui \prepi) |0)^ + sin /3j \prepj) \ 1)^ . (7) 

i j 

with \ai\' = = 1; \prepi), \prepj) denote tensor products of all qubit eigenstates 

but qubit z. 

The ground state (7) can be reached through quantum annealing, namely quantum 
ground state computationt^^'-*^^! . Since there is no output constraint, this amounts to solving 
a problem polynomial in network size. Therefore, for short, without entering into any detail, 
we assume that the time required to reach any desired probability of finding the network in 
the ground state (7) is polynomial in network size. 

In the right side of eq. (7), each term i satisfies the input constraints, all N.E., not the 
output constraint, whereas each term j satisfies also the output constraint and is therefore 
a network solution. 

sin^ 1? depends on the network and must decrease exponentially with network size, oth- 
erwise the problem would not be hard (solutions are "exponentially rare"). Of course: 

diagp, (0) = cos^ ^ |0), (0|, + sin^ ^ |1), , 

with sin^ -d "extremely small" . 

3) We operate now on qubit z as we did in Section IV, so that at some time r: 
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This is done by introducing the one-qubit Hamiltonian — E'^a^o^^O; with « E 
and ground state and by putting the qubit in interaction with a heat-bath whose 
temperature is lowered down close to zero. The internal energy of qubit z is thus (^^ if)) = 
(^f {t)\H^ I* (i)). At some time t^T ("after" relaxation), we have: {^^ (r)) ~ 0. 

As in Section IV, the network evolution can be split into two orthogonal components: 
\^{t))=a{t)\^' {t)) + P{t)\^" it))- 

• 1^' (t)) is ground (with reference to the network and the constrained input qubits 
internal energy, see Section IV), namely Ylih (0 1 -f^fe (^)) — 0) with z ranging 
over all the N.E. and the constrained input Hamiltonians; this state is hosted in the 
"good" Everett universe of amplitude a (t) ; 

(t) ) is excited and belongs to the subspace orthogonal to H\f , it is hosted in the 
"bad" Everett universe of amplitude /? (t) . 

We assume to be in the "good" Everett universe (until different notice). Qubit z relax- 
ation, under the continuous projection (i) and (ii) drives the transformation of the overall 

network state: 

I* (t)) = cos i^ + if (t)] ai \prepi) |0), + sin + if (t)] ^ f3j \prepj) ; 

« 3 

The driving condition is E^ cos^ [& + (p it)] = (^^ (t)); we should keep in mind that, in current 
assumptions, qubit r relaxation is independent of the rest of the network. Condition (ii) 
keeps all ccj and (5j unaltered throughout the evolution, as readily checked. At time t — t: 

j 

By measuring the network at time r, and by repeating the overall process for a sufficient 
number of times, a straightforward application of probability theory shows that the number 
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of repetitions required to check whether the network is satisfiable (with any desired prob- 
abihty that this check is correct) grows polynomially with network size. Furthermore, r is 
independent of network size. This would mean NP-complete = P. 

However, this result can be completely vanified if the probability of finding the network 
in the "good" universe decreases exponentially with network size. To find out how this 
probability decreases, we should compute how the qubit r-heat-bath interaction affects a {t) 
and (3 (t) . Doing this would require a level of formalization beyond the reach of the current 
study. 

Anyway, one ingredient of the "quantum computation speed-up" consists of substitut- 
ing interaction, which requires time, with interference, which requires no time. Particle 
statistics, in the current interpretation, works as interference and implements most of the 
network logical constraints. This, in a still generic way, should justify the hope of achieving 
a computation speed-up by using particle statistics. 

VII. DISCUSSION 

We have highlighted an alternative way of approaching quantum computation. With 
respect to time-sequential computation (Section I), it presents some significant differences: 

• logical constraints are simultaneous in time rather than being mapped on the time 
evolution of the computation process; 

• there is thus an analogy between such constraints and the simultaneous constraints 
(also logical in character) estabhshed by particle statistics; 

• in order to make this analogy explicit, we must introduce a blunt interpretation of 
particle statistics. A particle statistics symmetry would no more be a passive constant 
of motion which does nothing to a unitary evolution, but an active watchdog effect 
shaping that evolution; 

• the NP-complete SAT problem is in a way native in this approach. 
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The counterfactual interpretation of particle statistics might be interesting in itself. 
Checking it might not be out of reach. The link behaviour (Section IV) is similar to exchange 
interaction. One should devise an experimental situation, showing that an operation on one 
part drags the state of the other part according to the formulation propounded. 

We should finally note that this formulation of particle statistics can be represented in a 
two-way (advanced and retarded) propagation model. This has been done in ref. [12] and can 
be outlined as follows. We have been dealing with state vector reduction on a predetermined 
subspace, namely always on the constrained subspace 7i^"^\ never on the orthogonal subspace 
violating the constraint. Although elusively (since we are in counterfactual reasoning), this 
means an evolution affected by a condition coming from the future. 

Interestingly, a two-way propagation model of partial state vector reduction can also 
justify the speed-up of time-sequential computation. See ref. [12]. 

This research has been partly developed during the Elsag Bailey-ISI Workshop on Quan- 
tum Computation (Turin, 1997). Thanks are due to A. Ekert, D. Finkelstein, S. Lloyd and 
V. Vedral for useful suggestions. 
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